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GROUND-TRUTH  OBSERVATIONS  OF 
ICE-COVERED  NORTH  SLOPE  LAKES 
IMAGED  BY  RADAR 

W.F.  Weeks,  A.J.  Gow  and  R.J.  Schertler 


INTRODUCTION  radar  beam  actually  encounters  the  ice  at  angles 

that  are  appreciably  different  from  zero,  other 
In  a  recent  series  of  papers  (Sellmann  et  al.  factors  must  also  be  involved,  inasmuch  as  the 

1975,  Elachi  et  al.  1976,  Weeks  et  al.  1977, 1978)  bottom  of  undeformed  lake  ice  is  commonly  a 

attention  has  been  called  to  the  fact  that  the  very  flat,  specular  surface.  This  is  particularly 

radar  returns  from  ice-covered  lakes  with  appai-  true  for  the  thick,  undeformed  ice  covers  that 

ently  similar  upper  ice  surfaces  show  wide  varia-  develop  on  the  lakes  of  the  North  Slope  of 

tions.  The  patterns  of  these  variations  support  Alaska  Therefore,,  the  return  from  such  a 

the  hypothesis  that  the  areas  of  low  return  corre-  smooth  ice/water  surface  would  be  away  from 

spond  to  regions  where  the  ice  is  thick  enough  to  the  receiver  and  the  region  would  appear  as  a 

be  frozen  completely  to  the  lake  bottom,  result  low  return  in  the  SLAR  image,  even  if  a  high 

ing  in  an  ice/frozen  soil  interface.  Conversely,,  return  surface  was  present  (a  situation  similar  to 

where  a  high  return  is  observed  the  ice  cover  is  the  low  radar  return  observed  from  a  calm  sea), 

presumed  not  to  extend  completely  to  the  lake  Clearly,  for  a  strong  return  to  reach  the  receiver 

bottom,  resulting  in  an  ice/water  interface.  in  the  aircraft  additional  scattering  elements  are 

The  hypothesis  is  reasonable,  as  ice  and  fro-  required, 

zen  soil  have  nearly  the  same  dielectric  con-  It  has  been  suggested  (Weeks  et  al.  1978)  that 

stants  (a  low  contrast,  low  return  interface),  the  necessary  volume  scattering  within  the  ice  is 

while  ice  and  water  have  very  different  dielec-  the  result  of  the  numerous  highly  elongated  cy- 

tric  constants  (a  high  contrast,  high  return  inter-  lindrical  air  bubbles  that  commonly  occur  in  the 

face).  In  fact,  at  normal  incidence  the  reflection  ice  on  such  lakes.  Because  the  bubbles  are  elon- 

coefficient  at  X-band  frequencies  is  higher  by  a  gated  in  the  direction  of  heat  flow  (ice  growth), 

factor  of  about  1 2  for  an  ice/water  interface  than  the  elongation  in  an  undeformed  ice  cover  is  ex¬ 
it  is  for  an  ice/frozen  soil  interface.  The  differ-  actly  orthogonal  to  both  the  upper  and  lower  ice 

enc.es  between  the  reflection  coefficients  of  surfaces.  Accordingly,  the  bubbles  act  as  for- 

these  two  types  of  interfaces  could  possibly  be  ward  scatterers,  the  incoming  radar  pulse  is  only 

the  complete  explanation  if  the  radar  beam  was  scattered  downward  and  the  outgoing  radar 

perpendicular  to  the  ice/water  interface.  As  the 
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Figure  1.  Schematic  drawings  showing  the  nature  of  the  radar  scatter¬ 
ing  in  a)  lake  ice  with  elongated  bubbles,  b)  bubble-free  lake  ice,  and 
c)  lake  ice  that  is  frozen  to  the  bottom. 


pulse  (after  being  reflected  at  the  ice/water  inter¬ 
face)  is  only  scattered  upward.  Therefore,  for  a 
strong  radar  return  both  the  elongated  bubbles 
and  a  planar  ice/water  interface  are  required, 
while  a  weak  return  can  result  from  either  a  lack 
of  bubbles  or  an  ice/frozen  soil  interface.  If  the 
bubbles  in  the  ice  were  spherical  isotropic  scat¬ 
tered,  the  characteristics  of  the  lower  ice  inter¬ 
face  would  be  less  important  in  controlling  the 
strength  of  the  radar  return.  Figure  1  presents 
schematic  representations  of  these  different 
situations. 

Field  checks  of  these  hypotheses  are  limited. 
In  the  one  series  of  such  observations  that  has 
been  published  (Weeks  et  al.  1978)  six  sites  were 
visited  and  good  agreement  was  found  between 
the  observed  ice  conditions  and  the  nature  of 
the  radar  returns. 

In  the  present  paper  we  discuss  additional 
ground-truth  observations  made  at  a  number  of 
lake  ice  sites  on  the  North  Slope,  south  of  Lone¬ 
ly,,  Alaska,  and  compare  them  with  the  nature  of 
the  radar  returns.  We  also  utilize  a  new  set  of 
SLAR  images  to  prepare  maps  showing  areas  of 
the  commonly  shallow  North  Slope  lakes  where 
water  depths  are  greater  than  1.7  m  (the  approxi¬ 
mate  draft  of  the  ice  cover  at  the  time  the  radar 
images  were  obtained). 


THE  EXPERIMENT 

During  March  1979  the  NASA  Lewis  Research 
Center  C-131  aircraft  flew  a  series  of  missions  to 
study  the  radar  returns  from  sea  ice  in  the  Beau¬ 
fort  Sea.  On  the  way  to  and  from  the  sea  ice 
areas  it  was  possible  to  obtain  additional  radar 


imagery  of  a  portion  of  the  North  Slope  of 
Alaska  that  is  covered  by  a  myriad  of  shallow 
lakes.  As  part  of  this  area  was  included  in  some 
of  the  earlier  SLAR  studies  (Sellmann  et  al.  1975, 
Weeks  et  al.  1977, 1978)  comparisons  could  also 
be  made  with  previous  imagery.  The  general 
study  area  and  the  SLAR  coverage  are  shown  in 
Figure  2.  The  specific  SLAR  missions  utilized 
were  flown  on  16  and  20  March  1979.  The  qual¬ 
ity  of  the  SLAR  imagery  was  appreciably  better 
than  that  of  earlier  investigations.  The  character¬ 
istics  of  the  SLAR  system  used  are  given  in  Table 
1.  The  team  undertaking  the  ground-truth  obser¬ 
vations  (Cow  and  Weeks)  was  provided  with  high 
quality  SLAR  images  before  the  field  sampling 
was  undertaken.  In  retrospect  we  believe  that 
the  availability  of  imagery  during  the  field  obser¬ 
vations  was  essential  to  the  successful  selection 
of  suitable  field  test  sites. 

The  ground-truth  observations  were  made  on 
29  and  30  March  1979,  approximately  two  weeks 
after  the  SLAR  imagery  was  obtained.  As  the 
weather  during  the  intervening  period  had  been 
relatively  stable  (clear  and  cold),  there  was  no 
reason  to  expect  appreciable  differences  in  the 
nature  of  the  ice  cover  except  for  a  few  centime¬ 
ters  of  additional  ice  growth. 

Initially  we  planned  to  base  the  sampling  out 
of  Lonely,  traveling  back  and  forth  to  the  lakes 
with  a  snowmobile.  Fortunately,  at  the  last  min¬ 
ute  we  were  able  to  arrange  air  transportation  di¬ 
rectly  from  Barrow  to  the  sample  sites.  We  men¬ 
tion  this  as  we  are  now  of  the  opinion  that  it 
would  have  proved  most  difficult,  even  with 
good  maps  and  excellent  visibility,,  to  locate 
oneself  precisely  on  the  North  Slope  at  this  time 
of  year.  The  terrain  is  extremely  flat,  with  few 
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distinctive  features,  and,  in  addition,  the  entire 
surface  was  covered  with  snow.  During  the  two 
days  that  we  spent  collecting  samples,  visibility 
was  excellent,  and  low  sun  angles  helped  to  re¬ 
veal  slight  differences  in  surface  topography 
However,  we  were  commonly  forced  to  circle  a 
prospective  site  several  times  until  we  could 
identify  enough  features  to  fix  our  position. 
Once  we  had  landed,  accurate  positioning  usual¬ 
ly  was  impossible. 


Our  ground-truth  studies  of  the  lake  ice  fo 
cused  mainly  on  investigations  of  the  basic  phys 
ical  properties  of  the  ice.  These  studies,  made  or 
samples  obtained  with  a  CRREL  (7.6  cm)  coring 
auger,  included  measurements  of  density  elec¬ 
trolytic  conductivity,  and  bubble  and  crystal 
structure.  Coring  also  yielded  data  on  ice  sheet 
thicknesses  and  provided  access  to  the  water 
and/or  soil  underlying  the  lake  ice  covers. 


Table  1.  Characteristics  of  side-looking  airborne  radar  system. 

Type  Motorola  AN/APS-94D  Real  Aperture 

Frequency  X-Band,  Tunable  9  10-9.40  GHz 

Polarization  Horizontal-Horizontal 

Pulse  width  200  nanoseconds 

Output  power  45  kW  peak 

Pulse  repetition  750/second 

3-dB  antenna  beam  width  0.45  degree  =  7  85  mrad 

Range  resolution  30  meters 

Minimum  detectable  signal  -97  dBM 

Receiver  bandwidth  4  5  MHz 

Swath  width  50  km  each  side  of  aircraft 

Digital  data  processing 

Returns  from  transmitted  pulses  are  divided  into  3200  equal  right  antenna  and  3200 
equal  left  antenna  ground  range  bins.  For  each  range  bin,  the  SLAR  video  signal  is  sampled 
and  digitized  into  a  six-bit  data  word  The  digitized  SLAR  video  data  from  a  single  radar  im¬ 
pulse  are  accumulated  in  memory  and  averaged  exponentially  with  the  return  from  the  pre¬ 
vious  pulses  Auxiliary  data,  including  aircraft  drift,  ground  speed  and  location,  are  multi¬ 
plexed  with  the  average  radar  video  data 

Image  outputs  Real  time,  full  resolution,  analog  wet-processed  film 

Real  time,  digital  dry,  silver-process  paper  or  film 
Digital  tape  output  Instrumentation  tape  recorder  in  Bi-O-L  format, 

computer-compatible  tape  in  NRZ  format 


RESULTS  strongly  support  the  explanation  that  the  differ¬ 

ence  in  the  strength  of  the  return  is  associated 
Figure  3  is  a  map  showing  the  exact  locations  with  the  degree  of  dielectric  contrast  at  the  ice/ 

of  the  11  different  sites  sampled  in  1979  (sites  9-1  water  (high  contrast)  or  ice/frozen  soil  (low  con- 

to  9-11),  as  well  as  3  of  the  6  sites  sampled  in  trast)  interface.  Also,,  the  patterns  of  the  weak 

1976  (sites  6-A  to  6-C).  Figure  4  shows  the  radar  return  areas  correspond  well  with  areas  where 

imagery  of  the  same  area.  Table  2  gives  the  the  lakes  would  be  expected  to  be  shallow,  re- 

ground-truth  results  for  the  1979  sites  and  all  of  suiting  in  an  ice  cover  that  is  frozen  to  the  bot- 

the  1976  sites.  Also  indicated  is  the  nature  of  the  tom.  Examples  of  this  are  1)  the  dark  bands 

radar  return  at  the  site  in  question  (strong  or  around  the  edges  of  all  the  lakes  where  shallow 

weak).  The  ice  porosity  n  was  determined  ac-  water  has  to  occur,  2)  the  dark  band  separating 

cording  to  the  simple  relation  n  =  (e,-es)/et,  the  two  segments  of  Naluakruk  Lake  (Site  6-B), 

where  e,  and  e,  are  the  densities  of  pure  bubble-  which  is  the  former  land  between  the  two  lake 

free  ice  and  of  the  ice  sample,  respectively,,  segments  before  they  joined,  and  3)  the  long 

measi  ed  at  the  same  temperature.  The  electro-  dark  septum  extending  mto  Teshekpuk  Lake,  just 

lytic  conductivity  of  the  ice,  which  measures  the  to  the  west  of  Site  9-9,  which  also  represents  a 

bulk  content  of  dissolved  ions,  was  determined  submerged  boundary  between  two  former  lakes 

on  melted  samples  at  25°C,  while  the  conducti-  that  are  now  joined.  It  was  initially  intended  that 

vity  of  the  water  under  the  ice  was  measured  at  Site  9-9  be  located  on  the  septum.  However,  the 

ambient  conditions  (water  temperature  =  0°C).  aircraft  landed  too  far  to  the  east,  a  fact  verified 

An  examination  of  Table  2  and  Figures  3  and  4  on  take-off  when  both  the  sample  site  and  the 

shows  that  of  the  17  sites  sampled,  8  gave  weak  lake  outline  were  visible.  The  fact  that  the  water 

radar  returns  and  9  gave  strong  returns  All  sites  layer  under  the  ice  at  Site  9-9  was  only  0.34  m 

with  weak  returns  showed  no  water  under  the  ice  thick  is  also  interesting,  as  Teshekpuk  is  the 

(an  ice/soil  interface),  while  at  all  sites  with  a  deepest  lake  on  the  North  Slope  and  Site  9-9  is 

strong  return  an  ice/water  interface  was  present  well  away  (2.8  km)  from  the  nearest  land.  It  sup- 

There  were  no  exceptions.  These  observations  ports  the  interpretation  that  the  dark,  low-return 
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Figure  3.  Map  giving  the  exact  location  of  the  1 1  sites  sampled  in  1 979  (sites  9-1  to  9-1 1 )  and  three  of 
the  sites  sampled  in  1976  (sites  6 -A  to  6-Cj.  Also  indicated  is  the  nature  of  the  radar  returns  at  the  dif¬ 
ferent  lakes  (see  also  Figure  4).  Lined  lake  areas  are  regions  of  weak  returns  (depths  < 1.7  m);  unlined 
lake  areas  indicate  strong  returns  (depths  > 1.7  m).  Inspection  of  the  base  maps  used  to  prepare  this 
figure  (U.S.  Geological  Survey  Maps  Teshekpuk  C-1,  C-2,  D-1,  D-2  [1:63,360])  will  show  numerous 
small  lakes  and  ponds  not  indicated  in  the  figure.  These  all  showed  weak  radar  returns. 


septum  indicates  a  region  where  the  lake  is  even  bubbles  present  in  the  ice  during  the  late  winter 

shallower  and  the  ice  is  completely  frozen  to  the  of  1979?  Figure  5  presents  photographs  taken  in 

bottom.  transmitted  light  of  5-mm-thick  vertical  sections 

Tlie  other  component  of  the  present  explana-  of  ice  from  several  of  the  sampling  sites.  These 

tion  of  the  varied  radar  returns  from  lakes  with  core  sections  measured  7.6  cm  in  diameter.  Bub- 

apparently  similar  ice  (as  viewed  from  above)  is  bles  are  invariably  strongly  elongated  in  the  ver- 

that  the  elongated  bubbles  in  the  ice  give  the  tical.  As  is  obvious,  highly  elongated  (length/di¬ 
volume  scattering  necessary  tor  a  part  of  the  ameter  ratios  >25)  bubbles  are  common.  The 

radar  pulse  to  return  to  the  receiver.  Were  such  bubble  diameters  are  suprisingly  constant  at  a 
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Table  2.  Ground  truth  observations  from  lake  ice  sites  shown  in  Figure  a. 


Site  Nature  of 
no.  radar  return 


9-1  Strong 


9-2  Weak 


9-3  i 


9-5  Strong 


9-6  Weak 


Thickness  (m) 
Ice  Water 


1.88  0.77 


0.76 


•>93  0  37 


1  93  0  57 


i  67  0 


Density 

sample 

location 

(m) 

Ice 

density 

(Mg/m'J 

Ice 

porosity 

t%) 

Conductis'i > 
sample 
location 
(m) 

Conductivity 

(2SnC) 

(vmho/cm) 

1979  Sites 

013-0.21 

0.918 

0.0 

014-0  21 

12 

r  oo-i.io 

0.912 

06 

0  21-0.23 

4 

1.01-1  12 

10 

1  10-1 12 

16 

1  59-1.61 

15 

1  84-1.87 

24 

water 

600  (0°C) 

0  00-6  tt'J 

0  918 

0  0 

0  00-0.09 

13 

0«i-0.48 

0.885 

3.6 

0.09-0.11 

19 

-t» 

0  877 

4.5 

0  39-0.48 

25 

0  48-0.50 

16 

0  65-0.67 

335 

0  67-0.74 

540 

0.75-0  76 

847 

0  14-0,22 

0  902 

1.7 

015-0.22 

11 

0  26-C.44 

0.915 

0.3 

0.22-0  26 

19 

0  88-1.00 

C  317 

1.2 

0.26-0  34 

27 

0  34-0.43 

23 

0  44-0.48 

35 

0  88-0.99 

43 

1  04-1.08 

432 

.4i  0.96 

0.916 

02 

0.19-0  25 

12 

18-1.70 

0.904 

1  5 

0  26-0  28 

42 

0  82-0.84 

32 

0.84-0.89 

25 

0.89-0  95 

31 

0  96-0  98 

18 

1  23-1.25 

9 

1  58-1  69 

19 

1  70-1  72 

10 

water 

1350  (0°C) 

0.14-0.26 

0  918 

00 

014-0.26 

19 

1  07-1  20 

0  893 

2  7 

C.26-0  28 

5 

1  73-t  90 

0  902 

1.7 

1.05-1  07 

6 

1  07-1.20 

5 

1  73-1.75 

27 

1  90-1.94 

25 

water 

•500  (0°C) 

0  16-0.28 

0918 

00 

0.14-016 

20 

0  76-0.82 

0.889 

3.2 

0.17-0  22 

6 

1.06-1  15 

0.911 

08 

0.22-0  28 

22 

1.50-1.65 

0  915 

03 

0  82-0  84 

25 

1  06-1  15 

12 

1  15-1.18 

10 

1.50-1  57 

16 

1  57-1  64 

140 

1.65-1.67 

237 

7 


Table  2.  Ground  truth  observations  from  lake  ice  sites  shown  in  Figure  3  (cont'd). 


Density  Conductivity 


Site 

Nature  ot 

Thickness  (m) 

sample 

location 

Ice 

density 

Ice 

porosity 

sample 

location 

Conductivity 

I2S°C) 

no. 

radar  return 

Ice  Water 

(m) 

(%) 

(m) 

(fim  ho/cm) 

1979  Site*  (cont'd) 


9-7 

Strong 

1  92 

0.2S 

0,19-0  33 

0,918 

00 

0  19-0.26 

16 

0  83-0.96 

0.912 

0.6 

0  26-0.33 

29 

1.16-1.27 

0.889 

3.2 

0.33-0  36 

20 

1  68-1  83 

0  906 

1.3 

0.80-0  83 

44 

0  83-0.90 

8 

0  90-0.96 

30 

1.14-1.16 

9 

1  16-1.22 

7 

1  22-1.27 

32 

1.68-1  75 

91 

1.75-1.82 

37 

1.83-1.85 

62 

water 

1820  (0°C) 

9-8 

Weak 

1.70 

0 

water 

2600  (0°C) 

9-9 

Strong 

1  70 

0  34 

0.13-0  27 

0.915 

03 

0.11-0.13 

3 

0.98-1  10 

0.876 

46 

014-0.20 

14 

1.48-1  57 

0.898 

22 

0  20-0.27 

21 

0  96-0  98 

5 

0  99-1  05 

9 

1.05-1. TO 

30 

1  48-1  56 

12 

1  57-1.59 

7 

water 

30  (0°C) 

9-10 

Weak 

0.42 

0 

0  38-0  42 

0  916 

02 

0  00-0  03 

23 

0  03-0.05 

6 

0  27-0  35 

22 

0  35-0  38 

25 

9-11 

Weak 

0.15 

0 

0.00-013 

0  894 

2.6 

0  00-0.12 

1998 

013-0.15 

3740 

1976  Sites 

Water  salinity 
_  l°U 


6-A 

Strong 

2  29 

0 

0153 

6-B 

Weak 

2  21 

0 

6-C 

Strong 

2  31 

0  58 

0.163 

6-D 

Strong 

2.34 

0  66 

0  716 

6-E 

Weak 

0.71 

0 

6-F 

Strong 

2.29 

0.07 

0.406 

given  site.  Although  bubble  terminations  are  not 
rare,  they  do  not  appear  to  be  as  common  as  in 
lake  ice  from  more  temperate  regions.  That  this 
should  be  the  case  is  reasonable  as  the  ice  on 
North  Slope  lakes  undergoes  fewer  fluctuations 
in  growth  rate  during  its  formation  than  does 
lake  ice  in  more  temperate  regions.  Note  that 


the  porosity  of  the  ice  caused  by  the  bubbles 
(Table  2)  is  at  the  most  4.6%  and  is  usually  less 
than  3%.  Note  also  that  at  a  number  of  sites  the 
porosity  of  the  near-surface  ice  was  close  to 
zero.  These  measurements  are  supported  by  the 
field  observation  that  in  general  the  near-surface 
ice  contained  few  if  any  air  bubbles.  Apparently 
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an  appreciable  thickness  of  lake  ice  must  form 
before  the  water  below  the  lower  ice  surface  be¬ 
comes  sufficiently  supersaturated  with  gas  to 
cause  bubbles  to  nucleate.  We  suggest  that  the 
extremely  dark  areas  that  bound  the  edges  of  a 
number  of  the  low-return  lakes  (for  instance  the 
lake  containing  Site  9-2)  are  produced  by  near- 
bubble-free  lake  ice  freezing  to  the  bottom,  and 
that  the  slightly  higher  returns  from  the  center 
portions  of  such  lakes  are  caused  by  the  direct 
backscatter  of  the  radar  pulses  from  the  occa¬ 
sional  spherical  bubbles  that  occur,  as  well  as 
from  the  terminations  and  irregularities  in  the 
cylindrical  bubbles.  At  such  locations  all  the 
forward-scattered  radiation,  as  well  as  the  radia¬ 
tion  that  has  passed  through  the  ice  sheet 
without  encountering  a  scatterer,  will  penetrate 
into  the  soil,  resulting  in  a  low  return. 

The  impurity  level  of  the  majority  of  the  lake 
ice  is  low  (5  to  30  pmho/cm),  corresponding  to  ex- 
cdlent  quality  raw  water.  This  is  in  spite  of  the 
fact  that  the  impurity  levels  in  the  lake  water  are 
generally  high  (measured  conductivities  of  500 
to  1 820  pm  ho/cm  at  0°C,  corresponding  to  salini¬ 
ties  of  0  5  to  1  87„„).  (Contrast  this  with  a  typical 
New  Hampshire  lake,  with  water  conductivity 
between  40  and  60  pmho/cm  and  ice  conductiv¬ 
ity  frequently  less  than  2  pmho/cm  (Cow  and 
Langston  1977.)  The  reason  for  the  large  differ¬ 
ences  between  the  ice  and  water  compositions  is 
that  the  lake  ice  grows  with  a  planar  ice/water  in¬ 
terface,  and  the  amount  of  impurity  trapped  in 
the  ice  is  near  that  specified  by  the  equilibrium 
phase  relations.  The  exception  to  this  is  Site  9-11, 
where  the  water  being  frozen  was  sufficiently  sa¬ 
line  to  cause  a  non-planar  ice/water  interface  to 
develop,  resulting  in  the  entrapment  of  large 
amounts  of  impurities  in  the  ice.  The  ice  at  this 
site  was  sufficiently  saline  (2%o)  that  a  low  radar 
return  might  have  resulted  even  if  an  ice/water 
interface  had  been  present.  Another  tendency 
apparent  in  a  number  of  the  conductivity  pro¬ 
files  (Sites  9-1,  -3  and  -6)  is  the  pronounced  in¬ 
crease  in  conductivity  as  the  bottom  of  the  ice 
column  is  sampled.  This  is  especially  true  at  sites 
where  the  ice  is  frozen  to  the  lake  bed.  This  ef¬ 
fect  is  due  to  the  rapid  increase  in  the  impurity 
level  of  the  remaining  amount  of  residual  water 
as  the  result  of  the  rejection  of  impurities  by  the 
growing  ice. 

Other  than  at  Site  9-11,  where  the  ice  showed 
the  development  of  a  sea-ice-like  structure,  all 
the  ice  sampled  appeared  to  be  normal  congela¬ 
tion-type  lake  ice.  Structural  characteristics  of 
this  ice  are  best  demonstrated  in  thin  sections 


prepared  with  a  microtome  (a  histological  sec¬ 
tioning  device  modified  to  cut  ice)  and  photo¬ 
graphed  between  crossed  polaroids  to  reveal  the 
outlines  of  the  crystals  composing  the  ice.  Ice  at 
all  locations  was  first-year  ice,  that  is,  ice  pro¬ 
duced  during  a  single  winter. 

The  ice  texture  can  be  described  as  coarsely 
crystalline,  especially  in  the  deeper  parts  of  the 
ice  sheets  where  the  dimensions  of  crystals  in 
horizontally  sectioned  cores  frequently  exceed 
the  cross-sectional  diameter  of  the  cores  (7.6 
cm).  In  vertical  sections  of  core  the  texture  of 
the  ice  is  typically  columnar.  Orientations  of  c- 
axes  were  measured  on  sections  of  core  from 
several  lakes  in  order  to  obtain  some  idea  of 
overall  fabric.  Plots  of  c-axes  (Schmidt  net),  to¬ 
gether  with  the  crystal  texture  photographs  of 
ice  from  three  lakes,  are  presented  in  Figure  6. 
Though  most  thin  sections  contain  too  few  crys¬ 
tals  for  a  completely  adequate  determination  of 
the  fabric,  most  nevertheless  exhibit  a  horizontal 
distribution  of  c-axes— one  of  two  major  fabric 
patterns  observed  in  lake  ice  sheets.  Much  of  the 
ice  is  also  characterized  by  the  existence  of 
distinctive  subgrain  structure  in  many  of  the 
crystals.  Such  features  are  most  likely  produced 
by  stresses  induced  in  the  ice  lattice  during  the 
growth  of  the  crystals. 

At  Site  9-11,  where  the  lake  contains  a  high 
concentration  of  salts,  ice  with  a  structure  and 
fabric  closely  resembling  those  of  sea  ice  was 
formed  in  place  of  lake  ice.  A  thin  section  photo¬ 
graph  of  ice  from  near  the  ice/lake  bed  interface 
at  Site  9-11,  together  with  the  corresponding  c- 
axis  fabric,  is  presented  in  Figure  7. 

The  snow  cover  on  top  of  the  lakes  was  gener¬ 
ally  continuous,  although  a  few  scattered  patch¬ 
es  of  bare  ice  were  present  (Fig.  8).  Snow  depths 
were  usually  5  to  20  cm  and  occasional  small 
sastrugi  fields  were  observed.  Snow  densities 
were  estimated  to  be  in  the  0.3-  to  0.4-Mg/m3 
range.  There  was  no  evidence  of  any  correlation 
between  the  nature  of  the  snow  cover  or  its  dis¬ 
tribution  and  the  observed  differences  in  radar 
return. 

It  should  also  be  noted  that  there  are  irregular 
linear  regions  of  high  return  on  Naluakruk  and 
Teshekpuk  Lakes  and  to  a  lesser  extent  on  Okal- 
ik  Lake.  Similar  features  were  noted  but  at  slight¬ 
ly  different  locations  in  the  L-band  imagery  dis¬ 
cussed  by  Weeks  et  al.  (1978).  They  are  definitely 
pressure  ridges  caused  by  the  shifting  of  the  ice 
cover.  The  high  returns  are  attributable  to  the 
blocky  nature  of  the  ice  composing  the  ridges. 
Figure  9  shows  both  an  aerial  and  a  surface  view 
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Site  9-9 


2cm 


Figure  6  (cont'd). 
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Figure  7.  C -ax/s  orientation  plot  and  thir,  section  photograph  of  the  ice  from  Site 
9-11.  Note  the  sea-ice-like  texture. 


of  the  ridge  on  Teshekpuk  Lake.  At  places  the 
ridge  was  2  m  high. 

Finally,  one  perplexing  feature  in  the  imagery 
should  be  noted.  This  is  the  lower  return  ob¬ 
tained  from  the  western  portion  of  Teshekpuk 
Lake  The  dark  return  area  has  diffuse  bound¬ 
aries,  the  locations  of  which  change  with 
changes  in  the  direction  of  the  radar  overpass. 
Similar  but  less  pronounced  lower  return  areas 
are  also  observed  from  a  number  of  smaller 
lakes.  Such  changes  in  radar  return  do  not  ap¬ 
pear  to  be  related  to  the  bottom  topography  of 
the  lakes.  We  also  specifically  examined  the 
nature  of  the  snow  and  upper  ice  surface  to  see 
if  we  could  note  any  features  that  would  corre¬ 
late  with  the  low  return  area.  There  were  no  ob¬ 
vious  correlations. 


To  summarize  our  results  we  feel  that  our 
surface-truth  observations  have  provided  strong 
support  for  the  suggestion  that  the  primary  dif¬ 
ference  in  radar  return  noted  on  these  lakes  is 
the  result  of  the  presence  or  absence  of  water 
under  the  ice  with  an  ice/water  interface  giving  a 
high  return  and  an  ice/frozen  soil  interface  giv¬ 
ing  a  low  return.  The  field  observations  also  sup¬ 
port  the  suggestion  that  there  are  sufficient  elon¬ 
gated  bubbles  within  the  ice  to  provide  the  for¬ 
ward  scattering  required  if  there  is  to  be  a  return 
of  the  radar  pulse  to  the  receiver  on  the  aircraft. 
We  suggest  that  the  next  step  in  the  study  of  this 
general  problem  should  be  the  mathematical 
modeling  of  the  interactions  between  the  radar 
pulse,  the  upper  and  lower  ice  surfaces,  and  the 
elongated  bubbles  enclosed  in  the  ice. 
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Figure  9.  Aerial  and  surface  views  of  the  pressure  ridge  on  Teshekpu'  Lake. 


MAPS  OF  COMPLETELY  FROZEN 
NORTH  SLOPE  LAKES 

If  our  interpretation  of  the  radar  imagery  is 
correct,  then  existing  imagery  can  be  used  to 
prepare  maps  showing  the  lakes  or  parts  of  lakes 
that  are  frozen  completely  to  the  bottom.  We 
could,  of  course,  simply  reproduce  the  radar  im¬ 
ages  However,  they  would  be  very  large.  Also, 
the  image  degradation  and  scale  reduction  asso¬ 
ciated  with  preparing  such  illustrations  would  re¬ 
sult  in  the  loss  of  much  useful  detail.  Therefore, 
we  have  chosen  to  prepare  representative  maps 
that  summarize  the  nature  of  the  radar  return 
from  lakes  in  two  different  areas. 

The  average  measured  ice  thickness  at  the 
time  of  the  radar  overflights  was  1.88  m  (range 
1.70  to  1.93  m)  Assuming  that  the  ice  covers  are 
free-floating  and  of  approximately  constant 
thickness  throughout  the  study  area,  a  water 
depth  of  about  1.7  m  should  correspond  to  the 
boundary  between  regions  of  low  and  high  radar 
return  within  a  given  lake.  Figure  3  is  such  a  map, 
prepared  using  U  S  Geological  Survey  (1:63,360) 
topographic  sheets  Teshekpuk  C-1,  C-2,  D-1  and 
D-2  as  a  base  We  have  also  prepared  another  ex¬ 
ample  (Fig  10)  which  uses  the  Harrison  Bay  A-3 
topographic  sheet  as  a  base.  As  can  be  seen  in 
Figure  2,  the  area  covered  by  the  Harrison  Bay 
sheet  is  about  56  km  southeast  of  the  Teshekpuk 
Lake  area  The  blind  areas  in  Figure  3  are  too  far 
to  the  side  of  the  aircraft  to  provide  good  im¬ 


agery  (low  glancing  angles  and  poor  returns).  In 
some  locations  of  the  image  used  to  prepare  Fig¬ 
ure  10,  contrast  was  not  adequate  to  allow  a 
clear  determination  of  the  nature  of  the  return; 
such  lakes  are  marked  with  a  question  mark.  In 
addition,  the  imagery  in  the  vicinity  of  the  Col¬ 
ville  River  (Fig.  10,  lower  right  corner)  shows 
many  strong  returns  (presumably  due  to  the  pres¬ 
ence  of  steep  banks  and  gravel  bars)  which  make 
the  interpretation  of  the  returns  from  both  the 
river  ice  and  from  the  nearby  lakes  rather  spec¬ 
ulative.  Although  we  have  not  marked  this  on 
the  map,  most  of  the  Colville  gives  a  low  return, 
indicating  that  it  is  frozen  to  the  bottom.  A 
somewhat  different  problem  was  encountered  in 
the  upper  left  corner  of  Figure  10,  where  the 
steep  banks  and  narrowness  of  Judy  Creek  make 
interpretation  difficult. 

When  the  radar-return  maps  (Fig.  3  and  10)  are 
compared  with  the  original  USGS  topographic 
maps,  the  reader  will  note  the  actual  presence  of 
innumerable  small  lakes  that  are  not  indicated 
on  the  radar  maps.  All  of  these  small  lakes  give 
weak  radar  returns,  indicating  they  are  frozen  to 
the  bottom  In  using  Figures  3  and  10  it  should 
be  remembered  that  the  slant-range  radar  image 
contains  geometric  distortions  relative  to  the 
map.  In  preparing  the  map  we  have  attempted  to 
remove  these  distortions  visually.  Even  consider¬ 
ing  these  difficulties,  it  is  clear  that  the  SLAR  im¬ 
agery  and  the  maps  prepared  from  the  imagery 
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Figure  10.  Map  of  radar  returns  from  the  lakes  on  the  Harrison  Bay(A-3)  quadrangle 
(1:63,360).  Lined  lake  areas  are  regions  of  weak  returns  (depths  < 1.7  m);  unlined 
lake  areas  are  strong  returns  (depths  >1.7  m). 


contain  a  great  deal  of  useful  information  con¬ 
cerning  the  water  depths  of  North  Slope  lakes 
that  is  not  currently  available,  and  that  would  be 
expensive  to  obtain  by  conventional  means,  A 
person  interested  in  details  about  a  specific  lo¬ 
cation  is  advised  to  examine  the  original  imagery 
(available  at  CRREL  or  at  NASA  Lewis) 

In  conclusion  it  should  be  noted  that  in  gener¬ 
al  the  lakes  toward  the  Brooks  Range  (i.e.  toward 
the  south)  are  deeper  than  1.7  m,  while  the  lakes 
near  the  Arctic  Ocean  are  commonly  quite 
shallow.  The  imagery  also  indicates  that  .lakes 
giving  a  strong  return  are  particularly  rare  in  the 
vicinity  of  the  Prudhoe  Bay  oilfield.  If  this  had 
been  known  prior  to  the  development  of  the 
field,  it  would  have  been  helpful  in  site  planning 
as  well  as  in  prospecting  for  suitable  fresh  water 
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A  facsimile  catalog  card  in  Library  of  Congress  MARC 
format  is  reproduced  below. 
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